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Abstract The cluster statistics of BTW automata in the SOC states are 
obtained by extensive computer simulation. Various moments of the clusters 
are calculated and few results are compared with earlier available numerical 
estimates and exact results. Reasonably good agreement is observed. An 
extended statistical analysis has been made. 
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I. Introduction: 

Bak, Tang and Wiesenfeld (BTW) [1] have recently (1987) introduced 
the concept of Self Organised Criticality to describe the appearance of long 
range spatio-temporal correlations observed in extended, dissipative dynam- 
ical systems in nature. This phenomena of SOC is characterised by spon- 
taneous evolution into a steady state which shows long range temporal and 
spatial correlation. The simple lattice automata model of sandpile which 
shows this SOC behaviour was studied by BTW[1]. Substantial develop- 
ments have been made on the study of BTW model. Several properties of 
this critical state, e.g., entropy, height correlation, height probabilities etc 
have been calculated analytically. But the critical exponents have not been 
calculated analytically. Extensive numerical efforts have been performed to 
estimate various exponents. Recently the studies regarding this model have 
been reviewed by Dhar[2]. 

Here, in this paper, the statistics of the clusters, formed by the sites 
having particular value of the automaton, has been studied, by computer 
simulation. We have calculated the various statistical quantities e.g., average 
size of clusters, total numbers of clusters, maximum size of the clusters etc. in 
the SOC state and found the distribution function of the clusters. Related to 
the statistics of the clusters of the automation values, Manna has estimated 
the fraction of sites, having automation values, 0,1,2 and 3, present in the 
SOC state[3,4]. Majumder and Dhar calculated exactly the fraction of sites 
having minimum value of automation in the SOC state[5]. 

The paper is organised as follows: In section II the model and the simula- 
tion technique have been described, section III contains the numerical results 
obtained here, and the paper ends with a concluding section (section IV). 

II. The Model and Simulation: 

BTW lattice automata model [1] evolves to a stationary state in a self- 
organised (having no tunnable parameter) way. This state has no scale of 
length and time, hence called critical. Altogether the state is called self- 
organised critical state. The description of the lattice automata model is 
as follows : At each site of this lattice, a variable (automaton) z(i,j) is 
associated which can take positive integer values. Starting from the initial 
condition (at every site z(i,j) = 0), the value of z(i,j) is increased (so called 
addition of one 'sand' particle) at randomly chosen site of the lattice in 
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steps of unity as, 



z(i,j) = z(i,j) + 1. 



When the value of z at any site reaches a maximum z m , its value decreases 
by four units (i.e., it topples) and each of the four nearest neighbours gets 
one unit of z (maintaining local conservation) as follows: 



for z(i,j) > z m . Each boundary site is attatched with an additional site 
which acts as a sink (i.e., open boundary condition is used ). In this simula- 
tion, the value of z m = 4. 

It has been observed that, as the time goes on the average value (z) 
of z(i,j), over the space, increases and ultimately reaches a steady value 
(z c ) characterising the SOC state. At any time (i.e., a particular value of z), 
after the addition of one such 'sand' particle, when the system becomes stable 
(i.e., all the sites have z(i,j) < 4), the lattice contains the sites having the 
automaton values 0, 1, 2 and 3 only. These sites (having a particular value 
of z(i,j)) form clusters connected via the nearest neighbours. These clusters 
have some size distribution described by the function pk(sk), where pk(sk) 
denotes the number of clusters (formed via nearest neighbour connection of 
sites having their automaton value k) of size From this distribution Pk{sk), 
the various statistical moments of such clusters have been calculated here and 
studied as a function of z. The various statistical quantities (moments), for 
a particular class of cluster, are defined as follows: 

• Total fraction of sites (Fk(z)) having value k at a particular value of z. 



• Average size (Sk(z)) of clusters formed by the sites of automation value 
k. If Pk(sk) is the number of clusters (formed by the sites having their 
automation value z) of size s^. 



z(i,j) = z(i,j) -4 



z(i±l,j±l) = z(i±l,j±l) + l 



(1) 
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Total number of clusters (Nk(z)) formed by the sites having their au- 
tomaton value z(i,j) = k. 

N k{z) = J2pk( s k) 

Largest size (s™ ax {z)) of cluster formed by the sites of their automation 
value k. 

L k (z) = Max[s k ] 



III. Results: 

We have calculated the cluster size distribution, Pk(sk), in the SOC state 
for different automaton values (K = 0, 1,2 and 3). We did it by using the 
program [6] to calculate the cluster size distribution used in standard per- 
colation problem [7]. From the distribution all the statistical quantities (or 
moments) i.e., F k) S k , N k and L k are calculated in the SOC state. All these 
data has been obtained for L=400, by averaging over 20 different random 
samples. The values of F k have been estimated earlier [3,4,5] and our values 
agree well with the previous estimates (see Table I). 

Fig.l shows the semilog plot of the normalised cluster size distribution, 
(n s = p k (s k )/L 2 , is the number of clusters of size s per lattice site), for the 
automation value, k=3, and for L=400. The distribution function fits with 
the curve 

n s = .000004 * s -1 - 5 * exp(-.0A5 * s). 

To compare the cluster size distribution with the known result of perco- 
lation, we have plotted here the normalised cluster size distribution, n s = 
p k (s k )/L 2 , instead of p k (s k ). We see that this distribution function (for k 
= 3) has the same functional form as that of the percolation [7] (below p c ) 
apart from the exponent 9 of the pre-exponential factor, which may be the 
due to the fact that the sites are not randomly filled and they are correlated 
[9]. Another interesting point to note here is that the clusters (for k — 0, 
1, 2, 3) formed by automaton values are not system spanning (see the first 
column of table-I), although it can give a system spanning (scale invariant) 
avalanche in SOC state. It is to be mentioned that the other two distribution 
functions of the clusters, for automation value k=2 and 1, have also been ob- 
served. These two distribution functions are also of the same nature as that 
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of k=3. Since the maximum size of the cluster decreases as the automaton 
value decreases, the fluctuation in the distribution function for large clusters 
increases and hence these results have not shown explicitly. 

IV. Concluding Remarks: 

The statistical aspects of BTW cellular automata is investigated by com- 
puter silmulation. In the SOC state the automaton can take the values 0, 1, 
2 and 3. The cluster formaed by the fixed value of automaton (say 3) is iden- 
tified and measured its size. The clusters are defined via nearest neighbour 
connections. The statistical distribution of these clusters are calculated and 
from this distribution various moments are calculated in SOC state. Few 
results (for example the fraction of sites having a fixed value) are compared 
with the earlier numerical estimates (Ref.3-4) and exact calculations (Ref.5). 
The values of S^ oc and Lf° c for k=0 listed in table I are exactly known 
from Dhar's [2] burning algorithm. This study is a generalisation of earlier 
study (which estimates only fraction of sites having a fixed value) extending 
the calculation of entire statistics of BTW automata. It would be interest- 
ing to study the dynamical evolution (towards SOC state) of the cluster, its 
formation nucleation, coalescence etc. 
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Table I 



k 




qSOC 

k 




L soo 





0.0741 
0.0730 a 
0.0736 6 
0.0736 c 
0.07363 d 


1.000 


11856.2 


1.00 


1 


0.1746 
0.1740 a 
0.1740 6 
0.1739 d 


1.429 


19552.4 


10.85 


2 


0.3069 
0.307 a 
0.3062 6 
0.3063 d 


2.376 


20670.2 


37.90 


3 


0.444 
0.446 a 
0.4462 6 
0.4461 d 


4.326 


16430.5 


109.84 



a S. S. Manna, J. Stat. Phys., 59 (1990) 509 

b P. Grassberger & S. S. Manna J. de. Physique, 51 (1990) 1077 
C S. N. Majumder & D. Dhar, J. Phys. A:Math. Gen 24 (1990) 
L357 

d V. B. Priezzhev, J. Stat. Phys. 74 (1994) 955 
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Fig.l. The normalised distribution (n s ) of clusters of size s for the automa- 
ton value k — 3. Solid line represents the best fit. 
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